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Abstract 



This paper examines the parameter regimes in which coupled atomic and 
molecular Bose-Einstein condensates do not obey the Gross-Pitaevskii equa- 
tion. Stochastic field equations for coupled atomic and molecular condensates 
are derived using the functional positive-P representation. These equations 
describe the full quantum state of the coupled condensates and include the 
commonly used Gross-Pitaevskii equation as the noiseless limit. The model 
includes all interactions between the particles, background gas losses, two- 
body losses and the numerical simulations are performed in three dimensions. 
It is found that it is possible to differentiate the quantum and semiclassical 
behaviour when the particle density is sufficiently low and the coupling is 
sufficiently strong. 

PACS numbers: 03.75.Fi,03.75-b,32.80.Wr 
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I. INTRODUCTION 



One of the surprising discoveries that has been made since the experimental production 
of a Bose-Einstein condensate (BEC) in a weakly interacting gas [0] has been the fact that 
virtually all of their properties can be described using the Gross-Pitaevskii equation (GPE). 
This semiclassical approximation ignores any dynamics in the quantum statistics of the 
field. There has been much interest recently in the production of a molecular Bose-Einstein 
condensate (MBEC) from the photoassociation of an atomic BEC of a weakly interacting 
dilute gas We have recently calculated that in certain parameter regimes the GPE 

can give incorrect results for this process WM. This paper investigates the process with 
a three-dimensional model that includes background gas losses, two-body losses and the 
atom-molecule interactions. 



The success of the GPE |L5[ might not seem so surprising considering the effectiveness 
of the semiclassical approximation in quantum optics |T7J. It includes the effects of s-wave 
interactions, and can be readily generalised to include multi-component condensates with 
inter-species couplings [^|,|16| . As a semiclassical, mean-field theory it necessarily cannot give 
information about the quantum statistics of the condensates, but for most experiments with 
BEC these properties have not been observed. Quantum statistics affect some nonlinear 
quantum optical systems, the simplest of which is second harmonic generation, where pairs 
of photons are coupled to single, high-energy photons fL8| . The analogous process in atom 
optics is that of the coupling of a MBEC and a BEC, which may be done either through 



tuning of a Feshbach resonance ||11|| , or through photoassociation via a two-photon Raman 
coupling (§-|| . 

The Bose-enhancement of the photoassociation of atoms from a trapped BEC leads to 
giant, collective oscillations between the atomic and molecular populations. This enhance- 
ment of a chemical process was dubbed "superchemistry" by Heinzen et al. when they first 
modeled it using a two-component Gross-Pitaevskii equation (GPE) ||. A more recent 
model using the Hartree-Fock-Bogoliubov method includes pair correlations in the atomic 



field, and showed discrepancies with the results obtained with the GPE [11]. The full effects 



of the quantum nature of the fields can be modeled by a set of stochastic equations for 



the atomic and molecular fields based on the functional positive-P representation [19 22 



A recent paper used this technique in a one-dimensional calculation to show that it was 
possible to see the effects of the quantum statistics in the population dynamics of coupled 



atomic and molecular BEC WM - The signature was a reduction in the transfer from atoms 
to molecules. This paper uses a more complete model which includes all three dimensions 
as well as atom-molecule interactions and losses. The parameter regimes in which the GPE 
breaks down are determined. 

In Sec. H the extended model is described and the stochastic equations of motion for 
the system are derived. The following section is a brief examination of the behaviour of the 
GPE and the likely parameter regimes in which there are measurable deviations from this 
behaviour. In Sec. [TV] the evolution of the atomic population is compared with the solution 
of the GPE for a range of densities and coupling strengths. 
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II. MODEL 



An atomic field is coupled to a molecular field by two-colour Raman photoassociation. 
A first laser couples a single electronic level of the atomic field to a set of molecular excited 
states. A second laser then couples these states to a stable molecular level. The modes 
are arranged as shown in Fig. [I], with state |1) being the atomic BEC, state |2, v) the z^th 
vibrational level of the excited state of the MBEC and state |3) the stable MBEC Two 
laser fields induce a free-bound coupling between |1) and |2, v) and a bound-bound coupling 
between |2, v) and |3). In a rotating frame, the Hamiltonian may be written as 

i=l ij 

+ ihJ2 [ d 3 x (p v {x)ipl }l/ {x)ip 3 {x) - Ql{x)4) 2> „{x)i)l{x)) 
where ipj(x) is the field annihilation operator for the atomic or molecular field in state \j), 

A A th 

Tj and Vi are the kinetic and potential energy operators for the % field, Uij is the strength 
of the interatomic interactions between particles in states \i) and \j), k v (x) is the Rabi 
frequency of the free-bound photoassociation from level |1) to level |2,v) and Q(x) is the 
Rabi frequency of the bound-bound transition from level |2, v) to level |3). In this notation, 
the detunings of the lasers from the bare atomic and molecular energy levels are included 
in the potential energy terms Vi, v and V3. 
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FIGURES 



MBEC 



BEC 



FIG. 1. Energy level scheme for coherent free-bound-bound photoassociation. Levels |1), |2, v) 
and 1 3) are the electronic states for the atomic BEC, the z/-vibrational level excited MBEC and 
the stable MBEC respectively. 

In addition to the coherent effects produced by this Hamiltonian, the losses from the 
trapped levels are included by adding standard loss terms to the master equation. The total 
master equation is: 



P 



j 

+2 7 g } f d 3 xV[Mx)ij 3 (x)]p (1) 



where jj is the loss rate from level \j) due to background gases, jj is the two-body loss 

rate from level 7^3 is the two-body loss rate due to inelastic collisions between the BEC 
and MBEC fields, and where the superoperator T> is defined by 

V[c) = J[c] -A[c), (2a) 
J[c]p = cpc\ (2b) 

A[c]p = ^{c^cp + pc f c). (2c) 



A. Derivation of the stochastic field equations 

The master equation is a field operator equation with a non-trivial level of excitation, and 
is therefore impossible to solve numerically by direct means. Analytical results are precluded 
by the nonlinearities in the equations. To find a numerical method which is tractable 
in some parameter regimes, this master equation is written in the functional positive-P 
representation [|20|.p3|[ 

P({r,^},T) = p {a \{lft},r)\^^ ( 3 ) 
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where the density operator antinormally ordered with respect to the field operators in 

the Schrodinger picture. It is then possible to use the functional operator correspondences: 



■0p <-> ip a P, 



ftp 



[ft- 
ip a - 



d 



_d_\ 



P, 



(4) 



to write a functional Fokker-Planck equation (FPE) from the master equation. This FPE 
may be written in the form: 



dP 
~dt 



(5) 



where the elements fi and v correspond to the (4 + 2n) components of the fields in the 
positive-P representation: {ipf, tpi, . . . , 4>2 V , V'vj • • • > ^3 > V'f }> ^ i s the drift vector, D is the 
diffusion matrix, and d( v ++{ n ^y) = ^|y. The drift vector is given by 



A = 



(i)' 

2 



) + (-irn - if ] )ft^f + (-ir 13 - 7i 3 )^W? + E„ M>f^ v ^ 
-/d - ^) < + (ir n - 7 i 2) )^f^f + (iri 3 - 7i 3 )^^< + £„ «^ 



^3 - ^) ^ + (-ir 33 - i^W&f + Hr 13 - 713)^?^ - Ei/ W 2 a ,„ 



V 



^3 " ^ ) ^ + (^33 - 7?°) + (iFl 3 - 7l3K^V£ " E„ ^„ 



,32 



(6) 



where /Cj = —i/k(Tj + V}), and Tij = Uij/h. The terms proportional to the density of 
the excited molecular states have been dropped, as these upper states are going to be 
adiabatically eliminated. The diffusion matrix D is given by: 



/ -9n^f + E„ ^% 

-<?!*! 2 + E„ 



"#13^3^1 







"#1*3 



-gi 3 il>ZW \ 

-<7i* 3 VM 



"#33^; 



a2 o 



-9*33^ 



02 



(7) 



(2) 

where g, L j = 7^ + ir^- combines the atomic interactions and the two-body losses. The noise 
terms due to the self interaction of the excited molecular states have been ignored, as the 
density of these states is about to be assumed to be extremely small. This will enable them 
to be adiabatically eliminated. 
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The Fokker-Planck equation leads to the following set of equations for the excited molec- 
ular states: 

^ = -^f- fTi Al/ ^? - "T^V^ - y < + (8) 

= f-£v - + ^vVf ^ + ir 3i2 ,^3 a j < - y < + «:V* (9) 

In order for two lasers to provide a two-photon transition without populating the excited 
MBEC, the single photon detunings 5 V must be made very large. This means that the 
population of the excited level is very small and that its time derivative can be ignored. If 
the terms d t ip%„ an d Tip^u i n the equations of motion for ip^u are ignored on the basis that 
they are smaller than the other terms, it becomes possible to write solutions for ip^ixjt) 
explicitly in terms of "0? (x,t) and ip^ (x, t): 

= - 7$/ 2 - «Wfy? - ^ 3;2 ,^3 Q ) V£„ - y< + (10) 

o = (~v 2 ,, - 7S/2 + iTivA^i + •;<.•;;) <, - ^< + n* v ^. 

A rearrangement of these equations, and the further assumption that the single-photon 
detuning 5 V is larger than the interaction terms r^^iVf , the excited state loss rate 7^ 
and the trapping potential (i.e. V 2 , v ~ h5 u ), gives the following result: 

„ _ -i/2<^f 2 + z^ 
V>v - ^ • 

These equations can then be substituted into Eq. (||) and Eq. (0) to produce the drift and 
diffusion matrices for the two levels |1) and |3). It is then possible to find a matrix B such 
that D = B.B T . In general this matrix will not be unique. Finally, there is a theorem which 
allows us to map the solution of the master equation in Eq.(|I|) to the following set of Ito 
stochastic field equations: 

^ = Ui - ^ - (gx - «r)v>fy? - gi3 ^) ^ - ix^VS 

^ = ( -JCi - \- - (9l + «T)V>fy? - g{^t \ 4 + *X*^f 



+\/iX*Tpi-(gt+imi 2 V2 - vWf^C 2 (12) 
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where 



|2 



'v 



X = E^ (13b) 

A=5> !?y! (13c) 

and the 77^ are a set of real, Gaussian noise sources which are 5-correlated in time and space: 



^(x,0%(x',i') = 5 -5(x - x')5(t - 0, 
and are a set of complex, Gaussian noise sources which are (^-correlated in time and space: 



C(x, t)0(x', t>) = c%5(x - x')5(t - 0- 

These stochastic equations allow us to generate any multitime time-normally ordered 
quantum field averages by averaging selected moments of these fields over a sufficiently 
large sample of trajectories PU|. 



(rft(x, t)... &{x', t'yf&(x\ t") . . . v> V, r)) 

= 1pP(x, t) . . . 1pf3(x', t')lp a (x", t") . . . 1p a (x'", f"), 

where (T) is the (anti-) chronological time ordering operator. 

The Ito stochastic equations reduce to the GPE for this system if the noise terms are 
ignored. As can be seen later, this is not the case when using a different stochastic calculus, 
such as Stratonovich calculus, which reinterprets the noise terms and therefore requires 
corrections to the deterministic part of the equations. This reinforces the fact that the 
noise terms cannot be treated separately from the deterministic terms, as the same physics 
produces both of them. It is also worth reiterating that there is no useful interpretation of 
the individual trajectories of the simulation, only the ensemble averages. 

Our equations of motion are very similar to those obtainable by a direct coupling between 
the BEC and the stable MBEC The difference is the nonlinear light shift of the BEC, which 
is proportional to T, and the linear light shift of the MBEC, which is equal to A. 



III. SEMICLASSICAL SUPERCHEMISTRY 

Examination of the semiclassical equations of motion obtained by dropping the noise 
terms from Eq . (|l2|) shows that the behaviour is sensitive to the relative strengths of the 
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non-linearities. In the limit that there are no third order nonlinearities it can be seen that 
if the population starts in the atomic field then there is a complete, one-way conversion to 
molecules. In the presence of atomic and molecular interactions, the third order nonlinear- 
ities dephase this evolution, and cause a revival of the atomic population. This behaviour 
is shown in Fig. [| where the evolution is shown for different values of the atomic-molecular 
coupling rate. 




2 3 4 

time (ms) 

FIG. 2. Evolution of the semiclassical equations of motion for the atomic and molecular pop- 
ulations in a fixed three-dimensional harmonic trap. The values of x used were 3.0 x 10" 5 m 1 / 2 /s 
(solid), 3.0 x 1(T 6 m 1 / 2 / s (dashed), 3.0 x 10~ 7 m 1 / 2 /s (dash-dotted) and 3.0 x lO" 8 m 1 / 2 /s (dotted). 

This semiclassical behaviour has been observed in traveling-wave second harmonic gen- 
eration, which is a zero- dimensional analogue of coupled atomic and molecular condensates. 
However, it has also been shown that the semiclassical approximation gives incorrect predic- 
tions for the mean behaviour of the fields in traveling- wave second harmonic generation ||18|| , 
as well as for their quantum statistical properties [23J,2I|. This discrepancy is most pro- 
nounced when there is nearly complete conversion to the second harmonic, which occurs 
when the third order nonlinearities are very small [PEfl . It has since been shown that it is 
also possible to observe the effects of the quantum nature of the fields in a one-dimensional 
model of coupled atomic and molecular condensates \TA \. 

The breakdown of the semiclassical approximation has a simple explanation. The semi- 
classical equations do not include the spontaneous emission terms which allow a molecule 
to break into two atoms in the absence of an atomic field. To include these processes in 
the model the full stochastic equations of motion must be considered, which include all the 
physics of the quantum fields. The noise terms of the evolution mimic the dephasing effect 
of third-order nonlinearities, and cause a revival of the atomic population. 

For coupled atomic and molecular condensates in three dimensions, with losses due to 
background gas collisions and two-body collisions, it may be significantly more difficult to 
observe the quantum nature of the fields. In higher dimensions there is a larger volume 
of low density fields than in low dimensions, which in combination with the third order 
nonlinearities and the losses will tend to "blur out" the strong conversion to molecules 
which can be observed in the semiclassical theory for low dimensions. This effect can be 
counteracted by ensuring that the field is as dilute as possible, which makes the (second 
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order in the fields) coupling between the BEC and the MBEC relatively stronger than the 
(third order in the fields) interaction terms. In order to minimise the diffusion due to kinetic 
energy and further reduce the effect of the interactions, the Raman transition must be made 
as strong as possible. 

Making the Raman transition as strong as possible requires maximum laser intensities 
and a two-photon resonant coupling (A = A). The practical limits on the total strength of 
the coupling are the requirement that the excited molecular state must not have significant 
population, and the need to avoid inducing large, attractive atomic interactions through 
the non-linear light shift due to the laser-assisted photoionisation. Unless there is a low 
excited state population, spontaneous losses from that highly energetic state could disrupt 
the condensates, and the adiabatic approximation used to eliminate the upper level may 
break down. This limit imposes the condition that the single photon detuning 5 must be 
significantly larger than the Rabi frequencies of the individual laser beams. The second 
requirement means that the non-linear light shift should stay of the order of the repulsive 
interactions. If the Rabi frequency of the first laser is chosen such that T = Fu, then the 
nonlinear light shift will completely cancel the interactions between the atoms. There will 
still be repulsive interactions between the molecules and attraction between the atoms and 
the molecules. 

Dilute condensates can be created by evaporating below the BEC condition using strong 
straps, and then adiabatically lowering the trap strengths in all directions. The trap depth 
remains large compared to the chemical potential of the condensate even for extremely weak 



traps ||27|| . The field cannot be arbitrarily dilute, however, as the accurate measurement 
of the fields depends on a sufficient column density. Many current experiments use highly 
anisotropic "cigar-shaped" traps, which allow a high column density along the long axis. 
For a given aspect ratio uo x = uj y = Auo z and number of trapped atoms N, it is possible to 
calculate the transverse trap frequency required to obtain the required peak column density 

Or, 



' v 



18 al U\ x 7T 3 



Ux ~ 1,125 TV 3 m 2 A 2j ■ ( 14 ) 

The peak density scales inversely with the aspect ratio A, which has been as large as 100 
in current experiments . In the next section the behaviour of the quantum model of the 
system will be investigated in both dilute and strong traps, and with both weak and strong 
coupling. 



IV. QUANTUM SUPERCHEMISTRY 

This section will examine the difference between the quantum solution and the semi- 
classical solution allowing for physical constraints such as limited detuning and laser power. 
To ensure that the adiabatic approximation is being made self-consistently, the density of 
atoms in the excited state (which must remain small) can be calculated from the amplitudes 
of the other fields. 

Since our last calculation [H], the interaction between the atomic and molecular species 



has been measured by examining the resonances in the molecular formation due to pho- 
toionisation ||. The atoms are attracted to each other with a scattering length of 
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ai3 = —180 ± 150 a a , which leads to a value of U13 = 6 7ra 13 /?, 2 /m = —1.4 x 10~ 50 J m 3 . 
The strength of the interactions between the RD2 molecules has not been measured, so it 
will be assumed that they are the same as the interatomic interactions for the purposes of 
this calculation (Un = U 33 = 3 x 10~ 51 J m 3 ). 

The losses due to background gases are very slow on the timescale of the superchemistry. 
All calculations in this work use a value of 71 = = 0.01s -1 , which is a comfortable 



overestimate of the background gas losses measured by the group at JILA [29]. At the high 



densities involved in the production of an atomic BEC, the three-body losses are dominant 
and the two-body losses are negligible ('y^ < 1.6 x 10 -22 m 3 /s). By contrast, the three- 
body losses can be ignored at the low densities considered in this work, but in the presence 
of molecules the two-body losses cannot be ignored, as there is a significant cross-section 
for the inelastic scattering of an atom from a molecule. This was recently measured to be 
713^ < 8 x 10~ 17 m 3 /s, and in this paper the upper bound shall be used in all calculations. 

It is experimentally difficult to produce large Rabi frequencies for the atom-molecule 
interaction due to the low Franck-Condon factors, but this does not appear to be a limitation 
of this system provided large Rabi frequencies can be achieved for the molecule-molecule 
transition. 



A. Numerical methods 



Although Ito equations of motion for this system have been produced in the previous 
section, it is much easier to use high-order integration methods on Stratonovich equations. 
This is because the normal chain rule of differentiation applies to Stratonovich calculus, and 
standard high-order methods can be used without modification. The transformation from a 
set of Ito equations to a set of Stratonovich equations can be performed easily, as it involves 
a simple correction to the deterministic terms. Unfortunately, this correction is infinite for 
Eq. (|T2"D . This is a recurrence of the renormalisation problem often found in quantum field 
theories. 

The infinite correction only occurs in the continuous field description. When the con- 
version to Stratonovich calculus is done after the fields have been discretised and placed on 
a rectangular grid with spacing of Ax, Ay and Az in the three dimensions, the following 
correction terms are obtained: 



at 

\ <wL 1 

\ at ' 



( 



gi 3 /2+gi 
2 Ax AyAz 



1>l\ 
4 



correction 



2 Aa; AyAz 
_£l3/2+93_ / Q 
2 Ax AyAz ^3 
9*3/2+93 iJ3 
2 Ax Ay Az ^3 



(15) 



This correction is equivalent to a stepsize-dependent phase shift of the fields. It is clear 
that the correction is infinite in the limit of zero spatial cell size. In fact, this limit is not 
realistic for atomic fields, as the interactions between the particles are modeled as a contact 
potential, and this approximation fails at a sufficiently small distance scale. This puts a 
lower limit on the distance between the gridpoints, which will be of the order of magnitude 
of the scattering length. This forced discretisation is equivalent to the cutoff in momentum 
which is usually used to solve renormalisation problems. 
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In a numerical simulation with stepsize At, the noise terms ?y(x, t) are included at each 
timestep by choosing a random number R from a Gaussian distribution centred around zero 
and with unit width. The noise terms are then formed using r](x,t) = R/^Ax Ay Az At. 
This means that going to smaller grid spacing slows the computation in three ways. The field 
contains more points, which proportionally increases the computation time per timestep. 
Also, the noise terms become larger, which means that both a smaller timestep is required 
to perform a stable integration and a larger ensemble of paths is required to obtain precise 
averages. It follows that these stochastic methods are most successful when the density is 
low. Fortunately, the strong effect should also be in this limit. 

The stochastic integration was performed with the XMDS package, developed by Colle- 
cutt and Drummond f3(|. In collaboration with them and the ANU Supercomputer Facility, 
this has been adapted to run multiple integrations of the stochastic equations in parallel on 
the APAC supercomputer [EH] . 



B. Quantum superchemistry with varying coupling rates 

The feasibility of observing quantum field effects shall be demonstrated by considering 
the most dilute system possible. The diluteness will be limited by the requirement that the 
column density be large enough to provide easy measurements, and by the total number of 
atoms available in the initial trap. All further calculations assume that it will be possible 
to have a million atoms in the trap after the adiabatic expansion. The noise in the column 
density measurements is estimated to be 10 13 m -2 , which is that of the current BEC at 



ANU P2 |. For the purposes of this subsection, it is assumed that the trap has an aspect 
ratio of A = 30. The weakest trap in which 99.5% of the atoms are above the noise level in 
the column density has a trap frequency of uo x = 0.59s -1 , and a condensate size of 0.09mm 
in the strong axis and 0.3mm in the weak direction. The effect of varying the coupling rate 
in such a trap will be calculated in this subsection. The next subsection will examine the 
quantum field effects in stronger traps. 

Fig.(||) compares the evolution of the atomic population with the solution to the Gross- 
Pitaevskii equation. It can be seen that the effect of the quantum statistics of the field grows 
as the coupling increases. The coupling cannot be increased indefinitely, due to restrictions 
on both the available laser power and allowable detuning. 
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FIG. 3. Evolution of the atomic population (solid lines) compared to the solution of the GPE 
(dashed lines). Trap parameters are as given in the text. In each case the atomic- molecular coupling 
k was chosen such that T = Tn, and the detuning 5 was chosen to be 50 times larger than the 
molecule- molecule Rabi frequencies f2. Thus, a larger coupling strength requires a higher detuning. 
( a ) X = 4.7 x 10~ 6 m 3 / 2 s _1 , which corresponds to a detuning of 160MHz, (b) x = 1-5 x 10~ 5 m 3 / 2 
s -1 , which corresponds to a detuning of 1.6GHz, (c) x = 4-7 x 10~ 5 m 3 / 2 s -1 , which corresponds 
to a detuning of 16GHz. The dotted lines, where visible, are the uncertainties in the theory due 
to the sampling error. 

The difference between the dynamics of the system and the semiclassical result is essen- 
tially unaffected by the losses, which dominate at higher densities and over longer timescales. 
By contrast, the inclusion of the atom-molecule interactions in this model has reduced the 
distinction at high densities. From Fig.(|3]) it can be seen that the deviation from the semi- 
classical result is clearly observable for currently achievable trap and coupling parameters. 

C. Quantum superchemistry at varying densities 

Although it is possible to use very weak traps to work in the low density limit, most 
experiments are performed using stronger traps than those modeled in the previous section. 
This subsection will examine the effects of changing the density. 

The most dramatic effect of increasing the strength of the trap in the model is that 
the calculation becomes less stable. As previously explained, the noise in the numerical 
simulation becomes relatively larger as the unit cell size decreases. Stronger traps mean that 
the overall size of the condensate decreases, and so does the minimum allowable stepsize in 
the spatial grid. It therefore becomes harder to use the stochastic field equations to model the 
system for long times. The stochastic method will always produce results up to a certain 
point in time, but as the individual trajectories become unstable numerically it becomes 
impossible to integrate further. 

Under the constraint that there are just enough atoms to provide easily measurable 
results, the ratio of the second-order and third-order nonlinearities scales as l^ 3 , whereas 
the condensate volume will scale as u~ 2 . The difference between the full, quantum evolution 
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of the system and the semiclassical evolution depends on the ratio of the second-order and 
third-order nonlinearities, so based on this scaling the distinction may be measurable for 
trap strengths where the positive-P method cannot calculate the correct result. 

Fig. (|4]) shows the evolution of the atomic population for three different trap frequencies. 
The aspect ratio of the trap and the coupling rate are kept the same. The number of atoms 
in the trap is varied so that each trap contains the minimum number of atoms required to 
maintain 99.5% of them above a column density of a = 10 13 m~ 2 . An approximate method 
of calculating this is that the number of atoms N required to obtain a peak density of a p is 



N 



18 Ul 



a; 



1/3 



125 M 2 A 2 ^ 



(16) 



The figure shows a significant effect of the quantum statistics of the field even for much higher 
densities. Unfortunately, it also shows the difficulty in performing quantitative experiments 
at these densities, where the theoretical results are difficult or impossible to obtain by the 
stochastic methods used in this paper. In each plot, the integration cannot be performed 
far beyond the region shown. 
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FIG. 4. Evolution of the atomic population (solid lines) compared to the solution of the GPE 
(dashed lines). The parameters used are T = Tn, x = 4-7 x 10 -5 m 3 / 2 s _1 , which corresponds 
to a detuning of 16GHz. The trap has an aspect ratio of A = 30, and trap frequencies (a) 
u x = 2ir x 0.092Hz, (b) lo x = 2ir x 1.0Hz, and (c) tu x = 2ir x 10Hz. The dotted lines, where visible, 
are the uncertainties in the theory due to the sampling error. 

The importance of including the quantum statistical effects actually appears to increase 
with density. This is surprising, as we would expect the third-order nonlinearities to be 
the main source of dephasing of the conversion, and they are getting stronger with density. 
This may be due to the fact that the densities considered here are so low that the main 
dephasing is being caused by the kinetic energy term. If this is true, then there will be an 
optimal density which will maximise the parametric conversion rate compared to the kinetic 
energy terms, but minimise the two-body interaction terms with respect to the conversion 
rate. Detailed testing of this hypothesis is difficult due to the numerical difficulties, and is 
beyond the scope of this work. 
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V. CONCLUSIONS 



This paper has shown that it is possible to observe the quantum statistical effects of 
the field in coupled atomic and molecular BEC. The effect is most clearly visible in the 
low-density and high-coupling strength limits. 

Although it is a very successful model, the GPE cannot be applied to every system of 
coupled BECs. While it may seem reasonable to expect that the quantum statistics will 
tend to affect the multi-time correlations of the field rather than the mean field, our result 
shows it is sometimes also important to include them when considering the equations of 
motion for moments of the mean field. The signature of the breakdown of the GPE occurs 
in the simplest experimental observable - the total atomic and molecular populations. 

The calculations in this paper were numerically intensive. Although the positive-P rep- 
resentation contains the full description of the quantum field, it cannot be applied to every 
system of coupled BECs. Individually, trajectories do not have to behave in a physical 
fashion and they can become unstable over time in the absence of damping. The technique 
is therefore often only useful over short timescales, or for examining systems which are well 
described by the coherent state basis which underlies the description. This is not usually 
a natural basis of the atomic field, although it is often a good basis for the optical field. 
A possible counterexample may be a continuously pumped atom laser operating well above 
threshold, for which the positive-P representation may be tractable. It is also to be expected 
that the behaviour of such a device would depend critically on the quantum statistics. 
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